Local Pairing at [/-impurities in BCS Superconductors 
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We analyse here the role d-electrons on Anderson- U impurities play in superconductivity in a 
metal alloy. We find that phonon coupling at impurities counteracts the traditional effects which 
dominate T c suppression in the non-magnetic limit. In some cases, we find that non-magnetic 
impurities can enhance T c . Qualitative agreement is found between the predicted increase and the 
experimental data for VI- VI degenerate semiconductors doped with Tl or In. In the Kondo limit, 
a Fermi liquid analysis reveals that it is the enhancement in the density of states arising from the 
Kondo resonance that counteracts pair-weakening. 
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When a non-magnetic Anderson-C/ impurityB is placed 
in a superconductor, two distinct mechanisms can op- 
erate to suppress the superconducting transition tem- 
perature, T c . First, resonant scattering between the U- 
impurity and the conduction electrons leads to a broad- 
ening of the impurity levels. Such broadening increases 
the amplitude for binding conduction electrons on the 
impurity thereby inhibiting pair formationEI. Second, the 
on-site Coulomb repulsion leads to a weakening of the 
pairing interaction that keeps two electrocs. bound in a 
Cooper pair. As a result, T c is suppxesscclHu. 

In the non-magnetic limit, KondoEl impurities also lead 
to pair- weakening. When T < Tk, the formation of a 
Kondo singlet state at each impurity quenches the lo- 
cal momenta. However, the conduction electrons forming 
the many-body resonance around each impurity are spin- 
polarized. Consequently, conduction electrons of oppo- 
site spin experience a net Coulomb repulsion when they 
visit a Kondo impurityQ, thereby weakening the pair- 
interaction that holds a Cooper pair together. 

In .-.theoretical treatments of the pair-weakening 
eifectEr O, it is generally assumed that electrons on the 
impurities do not participate in superconductivity. That 
this view might not be entirely consistent can be seen 
from the early work of Ratto and Blandin (RB)tl Within 
an Anderson^/ model in a BCS superconductor, Ratto 
and Blandind showed that the Cooper pair amplitude 
on a U impurity is non-zero. Hence, electron pairs an- 
nihilated on a {/-impurity re-emerge in the conduction 
band as a Cooper pair. In addition, Suhl also suggested 
that local impurities should give rise to local regions of 
superconductivityE3. 

In this work, we consider explicitly phonon-induced 
pairing on non-magnetic Anderson-[/-impurities in a 
BCS superconductor. First, we show that the phonon 
coupling constants involving the impurity are at least 
as large as Xkk' , the standard phonon coupling constant 
for the Cooper pairs in the conduction band. As a re- 
sult, such local processes can Lead|-,to an enhancement, 
relative to previous treatmentsa'&B of T c . While it is 



well-known that pure japtential scattering can enhance 
T c in low-T c materialscij through coupling to transverse 
phonon modes, the present work suggests that in the case 
of non-magnetic [/-impurities, an additional channel is 
available to enhance T c . 

The starting point for our analysis is a collection of 
identical non-interacting (dilute limit) Anderson-i7 im- 
purities 
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In Eq. (1), ed is the defect energy of the impurity, Vik the 
overlap integral between a band state with momentum k 
and the i th impurity, a\ creates an electron in the band, 
a\ a creates an electron with spin a on the i th impurity 
and rii a — a irT ai a . In the Hartree-Fock limit, each impu- 
rity level is broadened with a width T = p (\Vik\ 2 ) where 
p is the density of states at the Fermi level. As a result 
of the hybridization of the localized level with electrons 
in the conduction band, the on-site Coulomb repulsion is 
felt by all electrons in the system. To include the pair- 
ing interactions in the superconducting state, we write 
the total Hamiltonian as H — H a + H pa i T where H pa i r 
contains the BCS interactions among all the electrons: 
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where the A's are determined by the electron-phonon in- 
teraction. The last two terms in Eq. (||) account for local 
pairing on the [/-impurity as well as scattering of Cooper 
pairs between the impurity and band states. In the non- 
magnetic limit, this problrm-i^as been solved previously 
without the last two terms 



1 



It is instructive at the outset to establish the magni- 
tude of the coupling constants in the last two terms in 
Eq. (0) . To evaluate Xd and Xik we expand the impurity 
states \io~) = on a ko\ko~) in terms of the k-states, \ka), 
in the band. In the expansion for the impurity states, 
we relied on the completeness of the k-basis. If the band- 
width D is finite, the k-states do not form a complete set. 
However, what is essential here is that the band contain 
the states |efc — €f\ < T. As is typically done, we assume 
that the matrix element Xkk' = A D is a constant for k- 
states with |ej. — ep\ less than ujp, the Debye frequency of 
the metal. In the estimates that follow, we will assume 
that u>d > r. Using the standard form for the electron- 
phonon interaction, V p h = X a' k+ *a\., _ q ,ak> ia,k-\ , we find 
that 



A d = <iT,U|V pfc |*T,U> 
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From the orthogonality of the k-states, it follows that linear combinations 
J2k \ a i<jka\ 2 = 1; hence, g(q = 0) = 1. From the con- 
tinuity of g(q) it follows that Xd = NX , where N is 
proportional to the number of electrons in the conduc- 
tion band. Hence, the on-site phonon interaction for the 
impurity electrons, is enhanced over the k-state pairing 
value. Consequently, the effective on-site Coulomb repul- 
sion is reduced to U = U + Xd- Similarly, the scattering 
matrix element 



in terms of matrix elements of the gap, 

At = — A D ^(afc'|a_fe/|) — £ ^{a^dii) 

k' i 

A., = -U{a^a,i{) - X ki ^(a fc / T a_ fe / x ) . (4) 



The Hartree-Fock on-site energy is E = €d + U(ni), with 
(rii) = (riii) = ( n il)- The presence of Xki causes the gap 
equations to become coupled. In fact, it is through this 
coupling that the single-particle density of states becomes 
enhanced. 

Let us define r\ = Xki / X a and introduce the Green func- 
tions G(p,q;t) = -{T[a ptT (t)al a (0)]) and F+(p,g;t) = 

(T alp^a^ (0) }. Here p or q represent either a local 

impurity or a band state. In terms of the discrete fre- 
quencies u> = uj n = (2n + l)7rT, the Fourier compo- 
nents of the Green functions are defined as G(p, q; t) — 
T XL z~ lu)t G u {p, q). The gap equations, Eq. (|) are then 
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is also related to A Q . An exact equality obtains if two con- 
ditions are true, namely (x\d) is real and Q-i^ki = &iiM- 
As we will see, the presence of the mixing term Xki en- 
hances the density of electron states participating in su- 
perconductivity. We will assume that Xki is a constant. 
Both effects, reduction of the on-site Coulomb repulsion 
and the enhancement of density of states at the Fermi 
level, play a positive role in superconductivity. We show 
ultimately that they can conspire to increase T c in the 
non-magnetic limit. 

A simple way to make these heuristic arguments rig- 
orous is through the Hartree-Fock decoupling of the 
Green function equations of motion method, vised by RB0. 
Whilcvpore sophisticated methods exist Jid'oEj the work 
of RBtl is sufficient to describe the non-magnetic limit of 
the Anderson model. The linearized Hartree-Fock equa- 
tions of motion for the creation operators can be written 
succinctly 



of the Fj| Green functions. The sum over k' in Eq.(^) is 
restricted over a momentum shell around the Fermi sur- 
face of width ujd- Eq. (||) must be solved to obtain T c . 
To facilitate this, we introduce the Hartree-Fock approx- 
imation to the Hamiltonian in the normal metal, H a as 
well as the corresponding Green function, G. From the 
Hartree-Fock equatio n of motion, (iuj — H )G tJ = 1 and 
the Gor'kov equationsllj, (iuj - -ff )G w + A Ft = 1 and 
(iui + Zf)F(t + A^G^ = 0, it follows that to linear order 
in the gap, F^(p,q) = G- u {-t, -p)A]<5 w (£, q), where t 
is summed over the k and i-states. This approximation 
is valid at and slightly below the critical temperature T c 
where the gap first appears. If we now substitute this 
expression into the self-consistent gap equations (Eq.(||)) 
and average over the random position of the impurities 
as well as average products of Green functions, we obtain 
a quadratic equation, 
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for the transition temperature where n s is the impu- 
rity concentration. We have introduced the average 
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Sw(p,q) = (G u (p,q)G- u (-p,-q)) av . In obtaining Eq. 
@ we decoupled the gap from the average of the prod- 
uct of Green functions. 

To facilitate a solution for T C) we note that the on- 
site repulsion U and the phonon coupling strength A Q 
are of quite different magnitudes. Typically, U 3> |A |. 
In this limit, Eq. (|^) simplifies to an equation linear in 
the phonon coupling, 
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where the subscript "eff" indicates division by 

The averages appearing in Eq. (0) can be evaluated 
straightforwardly following the ladder summation tech- 
niques. For example ,u 



^ Su(k, k ) ~ - — r tan - — - 



'E 2 + {\Lu\ + vy 



+ 0(nl) (8) 



The other averages are computed analogously. If we use 
these expressions for S u coupled with the standard BCS 
expression for the transition temperature, (|A |p ) _1 = 
In (2e7W£>/(7rT co )), we obtain that 
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where 



A = In (2WE 2 + r 2 /VT co ) - | tan" 1 | 
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1 + (U/ttE) tan- 1 (E/T) 
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with 7 is the Euler-Mascheroni constant. The corre- 
sponding expression for U c s can be obtained from Eq. 
( |l0| ) by replacing r\ with U. The local density on the im- 
purity, Pd{ep) =rT /( 7r (-£' 2 +r 2 ))j is given by the standard 
Lorentzian formEJ. Recall the 77 dependence arises from 
the scattering from a Cooper pair between the band and 
localized states. Also U < U as a result of the phonon 
coupling on [/-impurities. The importance of these terms 
should now be clear. When 7? e ff = and U — U, the cae.- 
rection to T c is precisely the negative correction of RB0. 
For rj ^ and U < U , the transition temperature is en- 
hanced relative-toJthe predictions in earlier treatments 
of this problemErEJ. In fact, we compare in Fig. (|l|) the 
predictions of the present theory for the initial slope of T c 
with the earlier predictions of RBcl. For modest values of 
r and U , we find that non-magnetic impurities can actu- 
ally enhance T c in contrast to the suppression indicative 



of pair-weakening. The magnitude of the increase in T c 
is of 0(n s eF / '(n T)) , where n is the conduction electron 
density. 
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FIG. 1. Theoretical values for the initial slope of T c pre- 
dicted from Eq. (^) as a function of the filling, (rid), on the 
impurity. p is the density of states, n s is the impurity con- 
centration and r = 2.0eV and U = 0.5eV. 

Experimentally, T c has been observed to increase-when 
transition metals were doped into TitZl. AndersonE£l has 
suggested that transition metals such as Fe are non- 
magnetic in Ti and hence might possibly increase T c . 
While the present theory is consistent with the exper- 
imental trends, the agreement should not be taken as 
a confirmation because the experimental samples con- 
tained unusually high dopant concentrations^. Further 
experiments are needed on such samples in the dilute im- 
purity regime to determine if non-magnetic impurities do 
in fact increase T c . However, in the context of degener- 
ate semiconductors such as PbTe and SnTe doped with 
Tl and In, respectively, the observed superconductivity 
has been attributed to arise solely from impurity statesEj. 
In SnTe doped with In, T c was increased by an order of 
magnitude with a 1% In-impurity level. More striking is 
the behavior in PbTe. In this material, superconductiv- 
ity with a transition temperature of T c — 1 — 2K was 
observed only upon doping with Tl. Dopants such as 
Na yield no superconductivity down to temperatures of 
T = .009K. Experimentally and theoreticallyE3, it is now 
well-accepted that local-phonon coupling at the dopant 
impurities is largely responsible for superconductivity in 
these semiconductors. In addition, the impurities are 
thought to be in the extreme mixed-valence regime as 
the on-sitejjepulsion is much less than the hybridiza- 
tion energyEJ. The large dielectric constant (e PbTe « 33) 
is primarily responsible for the lowering of the on-site 
Coulomb repulsion. For the experimentally relevant car- 
rier concentrations and an impurity doping level of 1%, 
we estimate that r^-ln ~ 1 and ep ~ .SeV. Also, T 
has been estimatedta to range between .01 to .leV. For 
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r = .leV, we estimate the magnitude of the relative 
increase in T c to be 0(n s eF/{n T)) as 10 which is quali- 
tatively consistent with the increase seen experimentally. 

We can extend this analysis to the non-magnetic limit, 
T < Tk, of the Kondo problem. In this limit (rii) = 1/2. 
Below Tk, a Kondo system is described by a a screened 
impurity in a Landau Feraii liquid with relatively weak 
quasi-particle interactions!]. Sakurau has shown that the 
non-magnetic limit of the Hartree-Fock treatment of an 
Anderson impurity can be used to describe a Kondo sys- 
tem for T < Tk by making the following transformation: 
1) r -► r/xn and 2) U cS -> = 7iTx Ti . We have 
introduced the vertex function Ti , for the inelastic scat- 
tering of a pair of d-electrons of opposite spin. Below Tk, 
the susceptibilities are given by \u = X|T = 7r r/(4T^). 
To calculate the transition temperature, we also need an 
expression for ri c g. According to Eq. (|9|), r\ and U are 
rescaled in the same way. Hence, in the ladder approxi- 
mation, the value of ?y c ff can be obtained by comparing 
two diagrams which correspond to (see Fig. (||a), and 
the diagram for the scattering of a pair of k-electrons into 
a pair of d-electrons as shown in Fig. (§b). We obtain 
that rjes = ri(iiT) 2 / (4TkU). Hence, the initial slope in 
T r is 



AT, 



A' - A' 



T co ^PoTk V U 
where A' = ln(8 7 (7V7r 2 T co )) - 1 
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FIG. 2. a) The vertex part of Y%\ in the ladder approxima- 
tion, b) the corresponding vertex for the scattering of a pair 
of k-electrons into a pair of d-electrons. 

To make contact with the Fermi-liquid picture of the 
Kondo problem, we rewrite this expression in the sugges- 
tive form 

AT r n. . I T K \ I ' . f Ti 



■In 



sen 2 " 



have dropped all irrelevant constants. Within the Fermi- 
liquid picture, T c cx Tij-exp(A ) where A is the dimen- 
sionless phonon coupling. In this expression, Tk replaces 
lod because electrons which are further away from the 
Fermi level than Tk are strongly scattered. In the pres- 
ence of [/-impurities, there are two corrections to the di- 
mensionless coupling constant A. First, we must include 
the enhancement in the density of states arising from the 
Kondo resonance. This enhancement!!] scales as tl s /Tk- 
In addition, we must include the repulsion between quasi- 
particle states of opposite spin. The repulsion-energy is 
essentially Tk below the Kondo temperatureQEil. Within 
the quasi-particle picture, this repulsion is spread over 
(p Tk) 2 states because there are two electrons partici- 
pating in each scattering event. Hence, the change in the 
dimensionless coupling constant is given byE3 
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However, AT C = —T C 5X/X 2 . Consequently, the initial 
slope in T c from the heuristic Fermi-liquid arguments 
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is identical in form to the more exact expression derived 
inEq. (O) because rj is 0(1). The second term in both of 
these expressions is the standard pair- weakening effect, 
whereas the first is a positive correction arising from the 
enhancement in the density of states at a Kondo impu- 
rity. In the strong-coupling regime, |A| > 1, Eq. (|l3| ) 
predicts that Kondo impurities can enhance T c . We con- 
clude then that non-magnetic impurities by virtue of lo- 
cal phonon pairing can counteract the standard T c sup- 
pressing effects and in some cases actually enhance T c . 
Experimental systems on which this prediction can be 
tested are the transition metal alloy Ti(Fc) and degener- 
ate semiconductors. 
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